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Abstract 

Expanding the ideas of the author's paper [25] we develop a pure game-theoretic 
approach to option pricing, by-passing stochastic modehng. Risk neutral probabili- 
ties emerge automatically from the robust control evaluation. This approach seems 
to be especially appealing for incomplete markets encompassing extensive, so to 
say untamed, randomness, when the coexistence of infinite number of risk neutral 
measures precludes one from unified pricing of derivative securities. Our method 
is robust enough to be able to accommodate various markets rules and settings in- 
cluding path dependent payoffs, American options and transaction costs. On the 
other hand, it leads to rather simple numerical algorithms. Continuous time limit 
is described by nonlinear and/or fractional Black-Scholes type equations. 
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1 Introduction 

Expanding the ideas of the author's papers [28], [29] we develop a pure game-theoretic 
approach to option pricing in a multi-dimensional market (rainbow options), where risk 
neutral probabilities emerge automatically from the robust control evaluation. The pro- 
cess of investment is considered as a zero-sum game of an investor with the Nature. 

For basic examples of complete markets, like binomial model or geometric Brown- 
ian motion, our approach yields the same results as the classical (by now) risk neutral 
evaluation developed by Cox-Ross-Rubinstein or Black-Scholes. However, for incomplete 
markets, like for rainbow options in mult i- dimensional binomial or interval models, the 
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coexistence of infinite number of risk neutral measures precludes one from unified pricing 
of derivative securities by usual methods. Several competing methods were proposed for 
pricing options under these circumstances (see e.g. a review in Bingham and Kiesel [T3]). 
most of them using certain subjective criteria, say a utility function for payoff or a certain 
risk measure. The difference in pricing arising from these methods is justified by referring 
vaguely to the intrinsic risk of incomplete markets. In our game-theoretic approach, no 
subjectivity enters the game. We define and calculate a hedge price, which is the minimal 
capital needed to meet the obligation for all performances of the markets, within the rules 
specified by the model (dominated hedging). 

Though our price satisfies the so called 'no strictly acceptable opportunities' (NSAO) 
condition suggested in Carr, Geman and Madan [Tl], one still may argue of course that 
this is not a completely fair price, as the absence of an exogenously specified initial 
probability distribution does not allow us to speak about a.s. performance and imphes 
per force a possibility of an additional surplus. To address this issue, we observe that 
together with the hedging price for buying a security, that may be called an upper price, 
one can equally reasonable define a lower price, which can be looked as a hedge for selling 
the security. The difference of these two values can be considered as a precise measure of 
the intrinsic risk that is incident to incomplete markets. An alternative way to deal with 
possible unpredictable surplus, as suggested e.g. in Lyons [39j for models with unknown 
volatility, consists in specifying a possible cash-back, which should be due to the holder 
of an option when the moves of the prices (unpredictable at the beginning) turn out to 
be favorable. 

Our method is robust enough to be able to accommodate various markets rules and 
settings including path dependent payoffs, American options, real options and transaction 
costs. Continuous time limit is described by nonlinear and/or fractional Black-Scholes 
type equations. 

As a possible weakness of our approach we should mention that, in order to be effective, 
eligible movements of a market should be reasonably bounded. Possible big jump should 
be taken into account separately, say by means of the theory of extreme values. 

Brief content of the paper is as follows. In Section |2] we set a stage by defining the 
game of an investor with the Nature leading to the basic game theoretic expression for 
the hedging price in the simplest case of a standard European (rainbow) option without 
transaction costs taken into account. 

In the next three sections, which are completely independent of any financial appli- 
cations, we carry out a preparatory work on evaluating certain rather abstract minmax 
expressions, showing in particular, how naturally risk neutral probabilities emerge (and 
more precisely the extreme points of these probabilities), as if by miracle, from minimizing 
Legendre transforms of concave functions defined on polyhedrons. 

In Section[6]we apply these results for the evaluation of hedge prices in the simplest set- 
ting. Section [7] shows the essential simplifications that become available for sub-modular 
payoffs. In particular, a unique risk -neutral selector can be specified sometimes, say in 
case of two colored options (for a still incomplete market). This is of crucial importance, 
as the major examples of real-life rainbow payoffs turn out to be sub-modular. Section 
|8] shows how transaction costs can be nicely fit into our model. Next two sections are 
devoted to the modifications needed for more complex models including path dependent 
payoffs, American options and transaction costs. Only in case of precisely J + 1 possi- 
ble jumps of a J-dimensional vector of stock process the corresponding market becomes 
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complete. 

Section [11] introduces the dual formulations and explicit expressions for upper and 
lower hedging prices. Next two sections are devoted to continuous time limits. These 
limits are obtained again without any probability, but only assuming that the magnitude 
of jumps per time r is of order r", a G [1/2, 1]. Finally, in Section HU the model with 
waiting times having power decay is discussed showing that its limit is described by a 
fractional (degenerate and/or nonlinear) version of Black-Scholes equation. 

Some bibliographical comments seem to be in order. Game-theoretic (or robust con- 
trol) approach for options was used in McEneaney [41], though in this paper the main 
point was in proving that the option prices of standard models can be characterized as 
viscosity solutions of the corresponding Hamilton- Jacobi equation. As a by-produce it 
was confirmed (similarly to analogous results in Avellaneda, Levy and Paras [3j and Lyons 
[39] ) that one can hedge prices in stochastic volatility models by the Black-Scholes strate- 
gies specified by the maximal volatility. A related paper is Olsder [16] , where only a basic 
one-dimensional model was analyzed, though with some transaction costs included. 

The reasonability of the extension of the binomial model allowing for price jumps inside 
the interval (interval model) was realized by several authors, see Kolokoltsov [28], Bernard 
[9], Aubin, Pujal and Saint-Pierre [2] and Roorda, Engwerda and Schumacher [19]. In the 
latter paper the term interval model was coined. The series of papers of P. Bernard et al [9] , 
[lOj . [TT] deals with one- dimensional models with very general strategies and transaction 
costs including both continuous and jump-type trading. Arising Hamilton- Jacobi-Bellman 
equation have peculiar degeneracies that require subtle techniques to handle. 

Hedging by domination (super- replication), rather than replication, is well establish 
in the literature, especially in connection with models incorporating transaction costs, see 
e.g. [1]. Problems with transaction costs in standard models are well known, as indicates 
the title 'There is no non trivial hedging portfolio for option pricing with transaction costs' 
of the paper Soner, Shreve and Cvitanic ^52j. This problem, similar to the story with 
incomplete markets, leads to the development of optimizations based on a subjectively 
chosen utility function, see e.g. Davis and Norman [15] or Barles and Soner [1]. 

Upper and lower values for prices were discussed in many places, see e.g. El Karoui 
and Quenez [23] or Roorda, Engwerda and Schumacher. [15]. An abstract definition of 
lower and upper prices can be given in the general game-theoretic approach to probability 
and finances advocated in monograph Shafer and Vovk [50j. 

The well known fact that the existing (whatever complicated) stochastic models are 
far from being precise reflections of the real dynamics of market prices leads naturally 
to the attempts to relax the assumed stochastic restrictions of models. For instance, 
Avellaneda, Levy and Paras [3] and Lyons [39j work with unknown volatilities leading 
to nonlinear Black-Scholes type equations (though still non-degenerate, unlike those ob- 
tained below). On the other hand, Hobson [23] (see also [22], [21] and references therein) 
suggests model independent estimates based on the observed prices of traded securities, 
the main technique being the Skorohod embedding problem (SEP). These approaches still 
build the theory on some basic underlying stochastic model (e. g. geometric Brownian 
motion), unlike our method that starts upfront with the robust control. Similarly, hedg- 
ing with respect to several (or all) equivalent martingale measures, based on the optional 
decomposition (see Follmer and Kramkov [TS] and Kramkov [36]), are based on some ini- 
tial probability law (with respect to which equivalence is considered). The risk- neutral or 
martingale measures that arise from our approach are not linked to any initial law. They 



3 



are not equivalent, but represent extreme points of risk- neutral measures on all possible 
realizations of a stock price process. 

'Fractional everything' becomes a popular topic in modern literature, see e.g. the re- 
cent monograph Tarasov [M] . For the study of financial markets, this is of course a natural 
step to move from the discussion of power laws in economics (see e.g. various perspectives 
in Uchaikin and Zolotarev [55j, Newman j45j, Maslov [40j and references therein) to the 
applicability of fractional dynamics in financial markets, see e.g. Meerschaert and Scala 
[13] . Meerschaert, Nane and Xiao [H], Jumarie [23], Wang [56] and references therein. 
Our game-theoretic analysis leads to degenerate and/or nonlinear versions of fractional 
Black-Scholes type equations. 

Notations. By |2;| we denote the magnitude (Euclidean norm) of a vector z and by 
11/11 the sup-norm of a function. We shall denote by o the point-wise multiplication of 
vectors (sometimes called i/adamard multiplication) : 

{y o zf = y'z\ 

Acknowledgements. The author is grateful to Sigurd Assing, Alain Bensoussan, 
David Hobson, Alex Mijatovic and Oleg Malafeyev for useful comments, and to Pierre 
Bernard for encouraging him to delve deeper in the topic of the paper. 

2 Colored options as a game against Nature 

Recall that a European option is a contract between two parties where one party has right 
to complete a transaction in the future (with previously agreed amount, date and price) 
if he/ she chooses, but is not obliged to do so. More precisely, consider a financial market 
dealing with several securities: the risk-free bonds (or bank account) and J common 
stocks, J = 1, 2.... In case J > 1, the corresponding options are called colored or rainbow 
options (J-colors option for a given J). Suppose the prices of the units of these securities, 
Bm and S*^, i G {1, 2, J}, change in discrete moments of time m = 1, 2, ... according to 
the recurrent equations B^+i = pB^, where the p > 1 is an interest rate which remains 
unchanged over time, and 5*^+1 = Cm+i^lny where G {1,2, ...,J}, are unknown 

sequences taking values in some fixed intervals Mi = [di,Ui] C R. This model generalizes 
the colored version of the classical CRR model in a natural way. In the latter a sequence 

is confined to take values only among two boundary points di,Ui, and it is supposed 
to be random with some given distribution. In our model any value in the interval [di, Ui] 
is allowed and no probabilistic assumptions are made. 

The type of an option is specified by a given premium function f of J variables. The 
following are the standard examples: 

option delivering the best of J risky assets and cash 

fiS\ S', S') = maxiS\ S', S' , K), (1) 
calls on the maximum of J risky assets 

f{S\ S\ S^) = max(0, max{S\ S\ S^) - K), (2) 
multiple-strike options 

fiS\ S', S') = max(0, S'-K,,S''- K2, S' - Kj), (3) 
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portfolio options 

f{S\ S^) = max(0, mS^ + na^' + ... + njS-^ - K), (4) 

and spread options 

f(S\ S^) = max(0, {S^ - S^) - K). (5) 

Here, the 5^, S*^, S*"^ represent the (in principle unknown at the start) expiration 
date values of the underlying assets, and K, Ki, Kj represent the (agreed from the 
beginning) strike prices. The presence of max in all these formulae reflects the basic 
assumption that the buyer is not obliged to exercise his/her right and would do it only in 
case of a positive gain. 

The investor is supposed to control the growth of his/her capital in the following way. 
Let Xm denote the capital of the investor at the time m = 1, 2, .... At each time m—l the 
investor determines his portfolio by choosing the numbers 7^ of common stocks of each 
kind to be held so that the structure of the capital is represented by the formula 

J J 

Xm-l = Yl <^Ll + (Xm-l - ^LSLl), 
j=l j=l 

where the expression in bracket corresponds to the part of his capital laid on the bank 
account. The control parameters 7^ can take all real values, i.e. short selling and bor- 
rowing are allowed. The value becomes known in the moment m and thus the capital 
at the moment m becomes 

J J 

if transaction costs are not taken into account. 

If n is the prescribed m,aturity date, then this procedures repeats n times starting from 
some initial capital X = Xq (selling price of an option) and at the end the investor is 
obliged to pay the premium / to the buyer. Thus the (final) income of the investor equals 

G{Xm S*^, S^, S^) — Xn — f{Sl^, S^, S^). (7) 

The evolution of the capital can thus be described by the n-step game of the investor 

with the Nature, the behavior of the latter being characterized by unknown parameters 
^Ij^. The strategy of the investor is by definition any sequences of vectors (71,- ■■ ,7„) 
such that each 7^ could be chosen using the whole previous information: the sequences 
Xq, X^_i and Sq, Sj^_i (for every stock j = 1, 2, J). The control parameters 7^ 
can take all real values, i.e. short selling and borrowing are allowed. A position of the 
game at any time m is characterized by J + 1 non-negative numbers Xjn, S^, • • • , S;^ with 
the final income specified by the function 

G{X,S\...,S')^X-f{S\...,S') (8) 

The main definition of the theory is as follows. A strategy 71, • • • , 7„, of the investor 
is called a hedge, if for any sequence (^1, • • • ,Cn) the investor is able to meet his/her 
obligations, i.e. 

G{Xn, S^, S^) > 0. 
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The minimal value of the capital Xq for which the hedge exists is called the hedging price 
H of an option. 

Looking for the guaranteed payoffs means looking for the worst case scenario (so 
called robust control approach), i.e. for the minimax strategies. Thus if the final income 
is specified by a function G, the guaranteed income of the investor in a one step game 
with the initial conditions X, S"^, S""^ is given by the Bellman operator 



. J J 

BGfX, 5\ ■ ■ ■ , 5^) = - max min G{pX ^S" -I'^S' - p^" -i'S\i'S\- ■ ■ S' 



J J 

. _ , . , _ >-^), 

(9) 

and (as it follows from the standard backward induction argument, see e.g. [S] or [53]) 
the guaranteed income of the investor in the n step game with the initial conditions 
Xo, 5*0, 5*0 is given by the formula 

B"G(Xo, 5*0, Sq). 
In our model G is given by ([H]). Clearly for G of the form 

G{X,S\--- ,S') = X-f{S\--- ,S'), 

and hence 

B«G(X, S\---,S')=X- \iB-f)iS\ ■■■,S'), 
where the reduced Bellman operator is defined as: 

{Bf){z\...,z') = mm max [/(^V, ^V, ■ ■ ■ , ^ V) - 7V(e - p)], (10) 

7 {(^e[d,,u,]} ^ 

or, in a more concise notations, 

i^f)iz) = min niax [f{^oz) - {-f,^oz-pz)]. (11) 

7 {V<^ldj,Uj]} 

This leads to the following result from [25] . 

Theorem 2.1. The minimal value of Xq for which the income of the investor is not 
negative (and which by definition is the hedge price if" in the n-step game) is given by 

H- = ^iB-f){Sl---,S^). (12) 

pn 

We shall develop a method for evaluating the operator ffTOj) . as well as its modifications 
for American options or when transaction costs are taken into account. 
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3 Underlying game-theoretic setting 



In this section we develop a general technique for the evaluation of minmax expressions of 
type ( ITOl) showing how naturally the extreme risk neutral probabilities arise in such evalu- 
ation. We also supply geometric estimations for these probabilities and the corresponding 
minimizing value of 7, which are crucial for a nonlinear extension given in Section [51 In 
order to explain the ideas clearly, we first develop the theory in dimension d = 2, and 
then extend it to arbitrary dimensions (which requires certain lengthy manipulation with 
multidimensional determinants). 

We shall denote by Int the interior of a closed set. Let a closed convex polygon in 
contains the origin as an interior point, and let ^1, ■ ■ ■ be its vertices, ordered 
anticlockwise. We shall denote such a polygon by 11 = n[,^i,--- ,^k]- The assumed 
condition 

OGlntn[ei,--- ,6] (13) 

implies that all do not vanish. 

We are interested in the following game-theoretic problem: find 

n[6,--- ,ed(/) = minmax[/(0 - (^,7)] (14) 
7GR2 Cen 

for a convex (possibly non strictly) function /. By convexity, this rewrites as 

mir--,W) = min max [/(e,) - (^,,7)]. (15) 

76R2 ii,---,(k 

Having this in mind, we shall analyze a slightly more general problem: for an arbitrary 
finite collection of non-vanishing vectors .^i, ■ ■ ■ ,C,k from R^, ordered anticlockwise, and 
arbitrary numbers /(^i), ■ ■ ■ , f{C,k), to calculate f lT^ (whenever the minimum exists). The 
corresponding polygon n[,^i, ■ ■ ■ ,C,k] (obtained by linking together all neighboring vectors 
^i, ^i+i, i = 1, - ■ ■ ,k, and ^k, ^1 with straight segments) may not be convex anymore. 

We shall start with the case of 11 being a triangle: 11 = n[^i,^2,^3]- Then condition 
( IT3|) implies that 7^ ^ct^j for a > and any Suppose the min in 

n[ei,6,e3](/) = min max [f{^,) - (^,,7)] (16) 

is attained on a vector 70 and the corresponding max on a certain ^j. Suppose this max 
is unique, so that 

/(e.)-(6,7)>m)-(0,7) (17) 

for all j 7^ i. As C,i 7^ 0, by changing 70 on a small amount we can reduce the l.h.s. of f lTTj) 
by preserving the inequality (fTTj) . This possibility contradicts the assumption that 70 is 
a minimal point. Hence, if 70 is a minimal point, the corresponding maximum must be 
attained on at least two vectors. Suppose it is attained on precisely two vectors, that is 

/te) - (e., 7) = f{Q - (0, 7) > fiu) - iu 7) (18) 

for some different i,j,m. Since the angle between is strictly less than vr, adding a 
vector 

e(e./ie,i + e,/ie.i) 
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to 7o will reduce simultaneously first two expressions from the l.h.s. of ( !T8|) . but preserve 
(for small enough e) the inequality on the r.h.s. of (fT8|) . This again contradicts the 
assumption that 70 is a minimal point. Hence, if 70 is a minimal point, it must satisfy 
the equation 

- (^1,7) = - (6,7) = /fe) - (^3,7), (19) 
which is equivalent to the system 

(6-ei,7o) = /(6)-/(ei), . . 

(e3-ei,7o) = /fe)-/(ei). ^ ^ 



Again by assumption (fT3|) . the vectors ^2 — ^1, ^3 — ^1 are independent. Hence system (120|) 
has a unique solution 79. 

For a pair of vectors m, f G R^, let f ) denote the oriented area of the parallelogram 
built on u, V and R{u) the result of the rotation of u on 90° anticlockwise. That is, for 

U = V = (f^f^), 

D{u,v) = u^v^ -u^v\ R{u) = {u^,-u^). 
Notice that the determinant of system fl20l) is 

^(6 - 6, 6 - 6) = ^(6, 6) + ^(6, ^1) + ^(^1, 6), 

and by the standard formulas of linear algebra, the unique solution 70 is 

^ /(eO^fe - ^3) + - + - 6) .21) 

D(6,e3) + /^(e3,ei) + ^(ei,6) ' ^ ^ 

and the corresponding optimal value 

nr. . _ /te)^(6,e3) + /(6)^(e3,6) + /fe)^(ei,6) ...^ 

Hence we arrive at the following. 

Proposition 3.1. Let a triangle H[^i,^2,^3] satisfy (IT^ . and let f {^1) , f {^2) , f {(,3) be 
arbitrary numbers. Then expression (1161) is given by (122 p and the minimum is attained 
on the single 70 given by fl2T]) . 



Proof. Our discussion above shows that if 70 is a minimum point, then it is unique and 
given by ( |2T]) . It remains to show that this 70 is in fact the minimal point. But this 
is straightforward, as any change in 70 would necessarily increase one of the expressions 
fi^i) ~ i^i^l) (which again follows from f lT3|) ). Alternatively, the same conclusion can be 
obtained indirectly from the observation that the minimum exists and is attained on some 
finite 7, because 

max [/(^,) - (^i,7)] 00, 
as 7 — )• 00. □ 
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Corollary 1. Expression (!22l) can be written equivalently as 

n[ei,6,e3](/) = E/(0, 

where the expectation is defined with respect to the probability law {pi,P2,P3} on ^i,^2,^3-' 



Pi 



D(6,e3) + i^(e3,6) + ^(6,6) 



((i,j,k) is either (1,2,3) or (2,3,1) or (3,1,2)). Moreover, this distribution is the unique 
probability on ^,1,^2,^3 such that 

3 

E{O = J2P.^^ = 0. (23) 

1=1 

Proof. Required uniqueness follows from the uniqueness of the expansion of ^3 with respect 
to the basis ^1, ^2- D 

We shall call a probability law on ^1,^2,^3 risk-neutral, if it satisfies ( 123|) . The reason 
for this terminology will be seen later. From the point of view of convex analysis this is 
just a probability on ^1,^2,^3 with barycenter in the origin. 

We can now calculate (IT^ for arbitrary k. 

Theorem 3.1. Let a polygon U = ■ ■ ■ ,^k] satisfy the following conditions: 
(i) No two vectors are linearly dependent; 

(a) The collection {^1, ■ ■ ■ ,C,k} does not belong to any half-space, i.e. there is no 
w G such that {u, ^i) > for all i. 
Then 

n[6, ■ ■ ■ ,W) = maxE,,^/(0 = me.x{p^f{^^)+P^fiQ+PTf(U), (24) 
where max is taken over all triples l<i<j<m<k such that 

OGlntn[e„ei,ed, (25) 

and {pf"^,pj"^,pl^} denotes the unique risk neutral probability on {^i,^j,^m} (given by 
Proposition \3.1]) with Ej^^ the corresponding expectation. 

Remark 1. Condition (i) is equivalent to the geometrical requirement that the origin does 
not lie on any diagonal ofU (or its extension), and condition (ii) is equivalent to (fT3|l . 



Proof. For any triple m} satisfying fl25|) . 

n[ei, ■ ■ ■ ,a](/) > min max [/(e) - (^,7)] = E^jmf{0, 

where Proposition 13.11 was used for the last equation. Hence 

n[ei,--- ,a](/) >maxE,,„,/(0. (26) 

A key geometrical observation is the following. Conditions (i) and (ii) imply that there 
exists a subset of the collection {^i, ■ ■ ■ ,^k} consisting only of three vectors {Ci,Cj,^m}, 
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but still satisfying these conditions (and hence the assumptions of Proposition 13. ip . This 
follows from the Caratheodory theorem (but can be also seen directly, as one can take 
an arbitrary and then choose, as Cm,^j, the vectors with the maximum angle (less 
than vr) with when rotating clockwise and anticlockwise respectively). This observation 
implies that the maximum on the r.h.s of (12^ is defined (the set of triples is not empty) 
and consequently the l.h.s. is bounded from below. Moreover, as for any triple {i,j,m} 
satisfying f l25|) . 

max [/(O - (^,7)] ^ 00, 

as 7 —7- 00, and hence also 

max [/(^i) - (^i,7)] 00, 
i=i,...,fc 

the minimum in (|T5l) is attained on some finite 7. Assuming that 70 is such a minimum 
point, we can now argue as above to conclude that 

f{^^) - {^^, 70) = /(O) - (0, 70) = fiU - iU 70) (27) 

for some triple ^i, ^j, ^m- Moreover, if these triple does not satisfy (i) and (ii), then (by the 
same argument) the l.h.s. of f l27|) can not strictly exceed f{^i) — (^z,7o) for all other ^i. 
Hence we are led to a conclusion that if 70 is a minimum point, then there exists a subset 
/ C {!,■■■ ,k} such that the expressions /(^i) — (6)7o) coincide for all / G / and the 
family I G /, satisfy conditions (i), (ii). But by the above geometrical observation, 
such a family has to contain a subfamily with three vectors only satisfying (i) and (ii). 
Consequently, (124|) holds. □ 

Remark 2. It is easy to see that the number of allowed triples {i,j,m} on the r.h.s. of 
(!24|) is two for k = 4, can be 3,4 or 5 (depending on the position of the origin inside U) 
for k = 5, and can be 4,6 or 8 for k = Q. This number seems to increase exponentially, 
as k ^ 00. 

Remark 3. Theorem \3.1\ can be easily extended to the situation when conditions (i) 
and/ or (ii) are not satisfied. Namely, if (ii) does not hold, then the l.h.s. of (fTSll is not 
defined (equals to —00). If (i) does not hold, then the max on the r.h.s of ([21]) should be 
over all eligible triples plus all risk neutral expectations over all pairs such that C,i = —<y^j, 
a> 0. 

Let us extend the results to higher dimensions d. Let us start with the simplest case 
of d + 1 vectors .^1, ■ ■ ■ , ^d+i in R"'- Suppose their convex hull n[^i, ■ ■ ■ , ^d+i] is such that 

OGlntn[ei,--- ,e<i+i]. (28) 
We are interested in evaluating the expression 

n[6, • ■ ■ ,ed+i](/) = min max[/(e,) - (^,,7)]. (29) 

A remarkable fact that we are going to reveal is linear in / and the minimizing 7 is unique 
and also depends linearly on /. 

Assume that R"^ is equipped with the standard basis ei, • ■ ■ , fixing the orientation. 
Without loss of generality we shall assume now that the vectors ^1, ■ ■ ■ ,^d+i are ordered 
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in such a way that the vectors {^2,^3, ■ ■ ■ form an oriented basis of R'^. The fact 

that the vector ^1 hes outside any half space containing this basis, allows one to identify 
the orientation of other subsets of ^1, ■ ■ ■ ,^d+i of size d. Namely, let {^j} denote the 
ordered subset of ^1, ■ ■ ■ , ^d+i obtained by taking out of it. The basis {^j} is oriented if 
and only if i is odd. For instance, if d = 3, the oriented bases form the triples {^2, ^3, ^4}, 
{6, 6, '^4}, {^i,'^4,6} and {^1,^3,6}- 

The same argument as for d = 2 leads us to the conclusion that a minimal point 70 
must satisfy the equation 



/(6)-(ei,7) = 
which is equivalent to the system 

te-ei,7o) = /te 



fi^d+i) - (6+1,7), 



d+l. 



(30) 



(31) 



From ( l28l) it follows that this system has a unique solution, say 70. 

To write it down explicitly, we shall use the natural extensions of the notations used 
above for = 2. For a collection of d vectors Mi, ■ ■ ■ ,Ud G R'^, let D{ui, - ■ ■ , Ud) denote 
the oriented volume of the parallelepiped built on ui,--- ,Ud and R{ui,--- ,Ud-i) the 
rotor of the family (mi, ■ ■ ■ ,Ud-i)- That is, denoting by upper scripts the coordinates of 
vectors, 



.Ud 



Ml 



ei det 



det 



Un 



u 



Un 



4/ 



R{ui, ■ ■ ■ ,Md_i) = det 



Ml 



ed\ 



-J 



u 



u 



62 det 



Ml 



Ml 



M 



d-l 



U 



+ 



d-l 



Finally, let us define a poly-linear operator R from an ordered collection 
of d vectors in R'^ to R*^: 



Ud] 



R{ui, ■■■ ,Ud) = R{u2 - Ml, M3 - Ml, 
= i?(M2, ■■■ ,Ud) - R{Ui,Us, ■■■ , Md) H h (-l)''"^i?(Mi, ■ ■ ■ , Md_i, 

Returning to system fl3T]) observe that its determinant, which we denote by D, equals 



Ud - Ml, 

\d-l 



D = D{^2-^i,--- ,^d+i-^i) = det 



\ ^d+l ^1 ^d+1 ^1 ■ ■ ■ ^d+1 ^1 / 

Using the linear dependence of a determinant on columns, this rewrites as 



^(6,--- ,ed+i)-ei'det 



d+l 



^2^ \ 



(^2 



-il det 



K^d 



d+l 



1 ^2^ 
Sd+1 



1 



it \ 



id+l, 
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implying that 



d+l 



(32) 



i=l 



Notice that according to the orientation specified above, D{{^i}) are positive (resp. 
negative) for odd i (resp. even i), implying that all terms in ( |32|) are positive, so that the 
collection of numbers 



i = 1, - ■ ■ ,d + 1, 



(33) 



define a probability law on the set ^i, ■ ■ ■ , C,d+i with a full support. 

By linear algebra, the unique solution 70 to system fl3Tl) is given by the formulas 



7o = ;^ det 



7o = ;^ det 



/(6)-/(6) e2-ef 

/(^d+i)-/(6) 

fe2-ei /(6)-/(6) 



and similar for other 7q. One sees by inspection that for any i 

, d+l 



d/ 



i=l 



and 



+ ■ ■ ■ + ;p (-i)''+'(/(ed+i) - /(^i)i?(6 - ^1, ■ ■ ■ , - ^i: 



which rewrites as 
1 



7o = - 



D 



fi^iMm) - /(6)^({6}) + ■ ■ ■ + {-irf{Ui)R{{Ui}). 

For example, in case c? = 3, we have 



7o 



-D234 + -D143 + -D124 + -D132 

/(ei)^234 + /(6)i^l43 + f{^3)Rl24 + f{U)Rl32 
D234 + -D143 + -D124 + -D132 



where Ajm = -D(Ci,0,^m) and 

-Rjjm = -R(6, ^j) + R{^j, ^m) + R{^m, 6)- 

As in case c? = 2, we arrive at the following. 



(34) 



(35) 



(36) 



70 = -^(/(6)-/(ei))^(e3-ei,--- ,ed+i-ei)-;^(/(e3)-/(ei))^(6-ei,e4-6,-- - ,ed+i-ei^ 



(37) 
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Proposition 3.2. Let a family {^i, ■ ■ ■ ,^d+i} in R*^ satisfy condition ( 128|) . and let /(^i), 
• • • , /(^d+i) &e arbitrary numbers. Then 

n[ei,--- = (38) 

i=l 

anc? i/ie minimum in fl29l) zs attained on the single 70 fl37|) . 

Corollary 2. Under the assumptions of Proposition l37^ 



n[ei,---,ed+i](/) = E/(0, 

m}) 



70 = E 



(39) 



(40) 



Dm)_ 

where the expectation is with respect to the probability law fl33l) . This law is the unique 
risk neutral probability law on {C,i, ■ ■ ■ ,^^+1}, i.e. the one satisfying 



d+l 



(41) 



i=l 



Proof. The only thing left to prove is that the law (l33l) satisfies ( HTl) . But as the r.h.s. of 
(HTl) can be written as the vector-valued determinant 



E/(0 = det 


/6 


6 • 

■ 


^d+1 






■ 


id+lJ 



that is a vector with co-ordinates 



det 



id+l/ 



. d. 



it clearly vanishes. 



□ 



To better visualize the above formulas, it is handy to delve a bit into their geometric 
meaning. Each term (— l)*^^D({i^i}) in fl32|) equals d\ times the volume of the pyramid 
(polyhedron) with vertices {0 U {^j}}. The determinant D, being the volume of the 
parallelepiped built on ^2 — ■Ci; ' ' ' ? ^d+i — (,1, equals d\ times the volume of the pyramid 
n[G) ■ ■ ■ ) ^d+i] in the affine space R'^ with vertices being the end points of the vectors ^j, 
i = I, - ■ ■ ,d + l. Consequently, formula (!32|) expresses the decomposition of the volume 
of the pyramid n[^i, ■ ■ ■ , ^d+i] into d+l parts, the volumes of the pyramids n[{0 U {^i}}] 
obtained by sectioning from the origin, and the weights of the distribution ( 133|) are the 
ratios of these parts to the whole volume. Furthermore, the magnitude of the rotor 
R{ui, ■ ■ ■ , Ud-i) is known to equal the volume of the parallelepiped built on ui, ■ ■ ■ , Ud-i- 
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Hence ||-R({6})II equals {d — 1)! times the volume (in the affine space R*^) of the {d — 1)- 
dimensional face of the pyramid n[^i, ■ ■ ■ ,^d+i] with vertices {^i}- Hence the magnitude 
of the ratios R{{^i})/ D{{C,i}), playing the roles of weights in (HOj) . are the ratios of the 
{d — 1)! times {d — l)-dimensional volumes of the bases of the pyramids n[{0 U {^i}}] to 
the d\ times their full d- dimensional volumes. Consequently, 

Dm) 

where hi is the length of the perpendicular from the origin to the affine hyperspace gen- 
erated by the end points of the vectors {^i}- These geometric considerations lead directly 
to the following estimates for expressions (!39|) and (HOj) . 



Corollary 3. 

|n[6,---,e.-fi](/)l< 11/11, (43) 

|7o|< 11/11. max /.ri, (44) 

1=1, ■■■ ,d+l 

with hi from (l42l) . 

These estimates are of importance for numerical calculations of 70 (yielding some kind 
of stability estimates with respect to the natural parameters). On the other hand, we 
shall need them for nonlinear extensions of Proposition 13.21 discussed later. 

Let us say that a finite family of non- vanishing vectors ^1, ■ ■ ■ , in R"^ are in general 
position, if the following conditions hold (extending naturally the corresponding conditions 
used in case d = 2): 

(i) No d vectors out of this family are linearly dependent, 

(n) The collection {^1, ■ ■ ■ ,^k} does not belong to any half-space, i.e. there is no 
u elV'' such that (w, ^i) > for all i. 



Remark 4. In Roorda, Schumacher and Engwerda JJBjj, condition (ii) is called positive 
completeness of the family {^1, ■ ■ ■ , ^k}- 

It is worth noting that in case k = d + 1, assuming (i) and (ii) is equivalent to (1281) . 
We are interested in evaluating the expression 

n[6, ■ ■ ■ ,a](/) = min max[/(e.) - (^.,7)]- (45) 

Theorem 3.2. Let a family of non-vanishing vectors ^1, ■ ■ ■ ,^fc in R"' satisfy (i) and (ii). 
Then 

n[ei,---,efc](/)=maxE,/(0 (46) 

where max is taken over all families I C {I,-- - , A;} of size \I\ = d + 1 that 

satisfy (ii) (i.e. such that the origin is contained in the interior o/H[^j, i G I]), and 
denotes the expectation with respect to the unique risk neutral probability on (given 
by Proposition \3.S\) . 

Proof. This is the same as the proof of Theorem 13.11 The key geometrical observation, 
that any subset of the family {^i,--- ,C,k} satisfying (i) and (ii) contains necessarily a 
subset with precisely d + 1 elements still satisfying (ii), is a direct consequence of the 
Caratheodory theorem. □ 
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Remark 5. As is easy to see, the max in (H6|) is attained on a family if and only 

f{i^)-{lI.i^)>f{^r)-{lI,ir) (47) 

for any i E I and any r, where 77 is the corresponding optimal value. Consequently, on the 
convex set of functions f satisfying inequalities ( H7|) for all r, the mapping ■ ■ ■ , C,k]{f) 
is linear: 

n[6,---,a](/) = E,/(0. 



4 Extreme points of risk-neutral laws 

We shall expand a bit on the geometrical interpretation of the above results. 

Let us call a probability law {pi, ■ ■ ■ ,pk} on a finite set {^1, ■ ■ ■ , ^k} of vectors in R'^ 
risk-neutral (with respect to the origin) if the origin is its barycenter, that is 

k 

E{O = Y,pS = 0. (48) 

i=l 

The geometrical interpretation we have in mind follows from the following simple 
observation. 

Proposition 4.1. For a family {C,i, ■ ■ ■ ,^k} satisfying (i) and (ii), the extreme points of 
the convex set of risk-neutral probabilities are risk-neutral probabilities with supports on 
subsets of size precisely d-\-l, satisfying themselves conditions (i) and (ii). 

Proof. It is clear that risk-neutral probabilities with supports on subsets of size precisely 
satisfying themselves conditions (i) and (ii) are extreme points. In fact, if this were 
not the case for such a probability law, then it could be presented as a convex combination 
of other risk-neutral laws. But these risk-neutral laws would necessarily have the same 
support as the initial law, which would contradict the uniqueness of the risk-neutral law 
supported on ci + 1 points in general position. 

Assume p = {p^,...,p"^) is a risk-neutral probability law on m > + 1 points 
^1) ■ ■ ■ ) ^m- Linear dependence of the vectors ^2 — ^i, - ■ ■ ,^m — (,i implies the existence of 
a non- vanishing vector b = (6^, . . . , 6™) in R"^ such that 

m m 

^6^ = 0, 5^6^^. = 0. 

i=l 1=1 

Hence for small enough e, the vectors p — eb and p + eb are risk neutral probability laws 
on ^1, ■ ■ ■ ,im- But ^ 

showing that p is not an extreme point. □ 

Proposition 14.11 allows one to reformulate Theorem 13.21 in the following way. 

Theorem 4.1. Let a family of non-vanishing vectors ii, - ■ ■ ,^k in R'^ satisfy (i) and (ii). 
Then the r.h.s. of formula fHB]) . i.e. 

n[6, ■ ■ ■ ,a](/) = min max[/(e,) - (^,,7)] = maxE,/(0 (49) 
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can be interpreted as the maximum of the averages of f with respect to all extreme points 
of the risk-neutral probabilities on ^i, ■ ■ ■ ,^k- All these extreme probabilities are expressed 
in a closed form, given by (!33|) . 

It is natural to ask, what happens if conditions (i) or (ii) do not hold. If (ii) does 
not hold, then n[^i, ■ ■ ■ ,C,k]{f) is not defined (equals to — oo). If only (i) does not hold, 
one just has to take into account possible additional extreme risk neutral probabilities 
coming from projections to subspaces. This leads to the following result obtained as a 
straightforward extension of Theorem 14.11 

Theorem 4.2. Let a family of non-vanishing vectors ^i, ■ ■ ■ R-*^ satisfy condition 

(ii). Then equation ( l49l) still holds, where the maximum is taken over the averages of f 
with respect to all extreme points of the risk-neutral probabilities on ^i, ■ ■ ■ ,^k- However, 
unlike the situation with condition (i) satisfied, these extreme risk neutral measures may 
have support not only on families of size d + 1 in general positions, but also on families 
of any size m + 1, m < d, such that they belong to a subspace of dimension m and form 
a set of general position in this subspace. 

Remark 6. Notice that min and max in fH9|) are not interchangeable, as clearly 

maxmin[/(^,) - ^,7)] = -00. 

Let us now formulate a mirror image of Theorem 14. H where min and max are reversed. 
Its proof is almost literally the same as the proof of Theorem 14.11 

Theorem 4.3. Under the assumptions of Theorem \5.1\ the expression 

n[ei,--- =max min [/(e,) + (e„ 7)] (50) 

can be evaluated by the formula 

n[ei,---,6](/)=minE,/(0, (51) 

where min is taken over all families / C {1, ■ ■ ■ , /c} of size \I\ = d+1 that satisfy 

(ii), and E/ denotes the expectation with respect to the unique risk neutral probability on 
The min in (ISTl) can be also interpreted as taken over all extreme points of the 
risk-neutral probabilities on ^i, - ■ ■ ,^k- 

Notice that the r.h.s. of ( I5T1) is similar to the formula for a coherent acceptability 
measure, see Artzner et al [1] and Roorda, Schumacher and Engwerda [48j. However, 
in the theory of acceptability measures, the collection of measures with respect to which 
the minimization is performed, is a subjectively specified. In our model, this collection 
is the collection of all extreme points that arises objectively as an evaluation tool for our 
game-theoretic problem. 

Remark 7. Coherent acceptability measures (j) introduced in Artzner et al |^ represent 
particular cases of nonlinear averages in the sense of Kolmogorov, see f35^ and Maslov 
14^ . The distinguished feature that leads to the representation of as an infimum over 
probability measures is its super- additivity. Clearly, postulating sub-additivity, instead of 
super- additivity, would lead similarly to the representation as a supremum over probability 
measures, and hence to the analog of ( H9|) . 
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5 A nonlinear extension 

Let us discuss a nonlinear extension of the above results. It will be used for the analysis of 
transaction costs. We start with the simplest additive perturbations, which are sufficient 
for the static (one-step) evaluations with transaction costs. 

For a finite set of non-vanishing vectors {^i, ■ ■ ■ ,C,k} in R'^, we shall evaluate the 
expression 

n[ei, • ■ ■ , ^k]if, g) = min max [/fc) - (6, 7) + 9il)], (52) 

where g is some continuous function. The main example to have in mind is 

g{x) = ci\x^\ H h Cd\x'^\ 

with some positive constants c,. We are going to make explicit the (intuitively clear) fact 
that if g is small enough, the min in (l52l) is attained on the same 7 as when g = 0. 

Theorem 5.1. Let {^1, ■ ■ ■ ,^k}, k > d, be a family of vectors in R'^, satisfying the general 
position conditions (i) and (ii) of Section\^ 

Let g{'y) be a non-negative Lipshitz continuous function that has well defined deriva- 
tives Dy[x) in all point x and in all directions y such that for any subfamily ^i, i & 
I C {^1, • ■ ■ ,ik}, which does not satisfy (ii), one can choose an u defining the subspace 
containing all ^i, i E I (i.e. (w,^i) > for all i E I) in such a way that 

{^i,oo)> D^{co), lel.^eR'. (53) 

Then the minimum in ( l52l) is finite, is attained on some 70 and 

n[ei,--- ,a](/,^7) = max[E,/(0 + ^7(7/)], (54) 

where max is taken over all families {^j}je/, / C {1, ■ ■ ■ ,k} of size \I\ = d + 1 that 
satisfy (ii) (i.e. such that the origin is contained in the interior ofI\.[^i, i G I]), and 
denotes the expectation with respect to the unique risk neutral probability on {^jjig/ (given 
by Proposition \3.S\) . and 'jj is the corresponding (unique) optimal values. 

In particular, if k = d + 1, then 70 is given by (!37|) . as in the case of vanishing g, and 

n[ei,--- ,ed+i](/,^7) = n[6,--- ,ed+i](/) + ^7(7o) (55) 

withn[iu--- ,ed+i](/) from dMD. 

Proof. Arguing now as in Section [31 suppose the min in (132]) is attained on a vector 70 
and the corresponding max is attained precisely on a subfamily ^j, i e / C {^1, ■ ■ ■ 
so that 

f{i^)-{lO,^i) 

coincide for alH G / and 

/(6) - (7, 6) + 9{l) > fi^j) - (7, ^j) + gil) (56) 

for j ^ I and 7 = 70, but this family does not satisfy (ii). (This is of course always the 
case for the subfamilies of the size |J| < d -\- 1.) Let us pick up an u satisfying (!53|) . As 
for 7 = 7o + ew, 

/(e.) - (e., 7) + gil) = fi^^) - 7o) + g{lo) - tu) - D^gi^o)] + o(e). 



17 



this expression is less than 

/te)-te,7o)+^7(7o) 

for small enough e > and all i E I. But at the same time ( |56l) is preserved for small e 
contradicting the minimality of 70. Hence, if 70 is a minimal point, the corresponding max 
must be attained on a family satisfying (ii). But any such family contains a subfamily 
with d+1 elements only (by the Caratheodory theorem). 

Do go further, let us assume first that k = d + 1. Then a possible value of 70 is unique. 
Moreover, the minimum exists and is attained on some finite 7, because 

max [/(e.) - (e., 7) + 9{l)] ^ 00, (57) 

as 7 — 00 (as this holds already for vanishing g). And consequently it is attained on the 
single possible candidate 70. 

Let now k > d + Ihe arbitrary. Using the case = d + 1 we can conclude that 

n[ei, ■ ■ ■ , ik]{f. 9) > max[E,/(0 + g{ii)l 

and hence the l.h.s is bounded below and fl57|) holds. Hence the minimum in fl52|) is 
attained on some 70, which implies flM|) due to the characterization of optimal 7 given 
above. □ 

Remark 8. In case d = 2, k = 3, condition (153!) is fulfilled if for any i,j and 7 G 

2\D,{^,)\ < |e.|max(|e.|,|Ol)(l + cos0(6,O)), (58) 

where by (j){x,y) we denote the angle between vectors x,y. 

Let us turn to the fully nonlinear (in 7) extension of our game-theoretic problem: to 
evaluate the minmax expression 

n[6, ■ ■ ■ = min max [f{^,,^) - (6,7)]- (59) 

Let us introduce two characteristics of a system {^1, ■ ■ ■ ,^k}, satisfying the general 
position conditions (i) and (ii) of Section El that measure numerically a spread of the 
elements of this system around the origin. 

Let Xi = >ii{C,i, - ■ ■ ,^k) be the minimum among the numbers x such that for any 
subfamily ^j, z G / C {^1, ■ ■ ■ ,^k}, which does not satisfy (ii), one can choose a vector 
Ui G of unit norm such that 

(ei,a;)>x, I el. (60) 

This xi is clearly positive by conditions (i), (ii). Let X2 = X2(^i, ■ ■ ■ , ^fc) be the minimum 
of the lengths of all perpendiculars from the origin to the affine hyper-subspaces generated 
by the end points of any subfamily containing d vectors. 

Theorem 5.2. Let {^1, ■ ■ ■ k > d, be a family of vectors in R'^, satisfying the general 
position conditions (i) and (ii) of Section \B 

Let the function f be bounded below and Lipshitz continuous in 7, i.e. 

l/te,7i)-/te,72)| <x|7i-72| (61) 
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for all i, with a Lipshitz constant that is less than both k\ and K2: 



X < min(xi, X2). (62) 
Then the minimum in ( |59l) is finite, is attained on some 70 and 



n[ei,--- ,ed(/) = max[E,/(e,7/)], (63) 

•where max is taken over all families I C {1, ■ ■ ■ ,k} of size \I\ = d + 1 that 

satisfy (ii), E/ denotes the expectation with respect to the unique risk neutral probability 
on {^i}i<zi (given by Proposition 1 3. and 7/ is the corresponding (unique) optimal value, 
constructed below. 

In particular, if k = d + 1, then 70 is the unique solution of equation fl65l) below. 

Proof. As in the proof of Theorem 15 -H using now (16T!) and (!60|) . we show that the mini- 
mum cannot be attained on a 7 such that the corresponding maximum is attained only on 
a subfamily i E I C ■ ■ ■ ,C,k}, that does not satisfy (ii). And again we conclude that 
if 7o is a minimal point, the corresponding max must be attained on a family satisfying 
(ii) and containing d + 1 elements only. 

Let us assume that k = d + 1. Then a possible value of 70 satisfies the system 

(e.-6,7o) = /(e.,7o)-/(ei,7o), i = 2,---,d+l, (64) 
which by (140|) rewrites as 



70 = E 



(65) 



where the expectation is with respect to the probability law fl33|) . This is a fixed point 
equation. Condition fl6T]) . fl62|) . the definition of X2 and estimate (jH]) imply that the 
mapping on the r.h.s. is a contraction, and hence equation ([H5D has a unique solution 70. 
Moreover, the minimum in fl5I?]) exists and is attained on some finite 7, because 

max [fi^i, 7) - {^i, 7)] -> 00, (66) 

as 7 — )■ 00 (as this holds already for vanishing /). And consequently it is attained on the 
single possible candidate 70. 

Let now > ci + 1 be arbitrary. Using the case A; = ci + 1 we can conclude that 

n[ei,---,a](/)>maxE,/(e,7/), 

and hence the l.h.s is bounded below and fl66|) holds. Hence the minimum in fl59|) is 
attained on some 70, which implies (|63|1 due to the characterization of optimal 7 given 
above. □ 

In applications to options we need to use Theorem 15.21 recursively under expanding 
systems of vectors ^. To this end, we require some estimates indicating the change of 
basic coefficients of spread under linear scaling of all co-ordinates. 

For a vector z G R'^ with positive coordinates let 

6{z) = max 2;*/ minz*. 
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Proposition 5.1. Let a system {^i, ■ ■ ■ ,^a:} of vectors in R'^ satisfy the general posi- 
tion conditions (i) and (ii) of Section Let xi,X2 be the characteristics of the system 
{^1) ■ ■ ■ introduced above and, for a vector z G R'' with positive co-ordinates, let 
>Ci{z), >i2{z) denote the characteristics of the system {z o ^i^ - ■ ■ , 2; o ^^}. Then 

Mz) > \z\Hi{d5{z))-\ X2(z) > \z\H2iVd5{z))-\ (67) 

Proof. Let us denote by z~^, just for this proof, the vector in R'^ with co-ordinates z^^. 
For a unit vector = |2;~^|~^2~^ o co, we get, using > k\ that 

{zo^i^^) = \z~'\{ii,Uj)>\z\Kx 



z\ \z~^ 



Hence to get the first inequahty in fl67|) it remains to observe that 

-'\ = {Y^z]Y^zffi^<d^{z). 



\z\ \z 



Turning to the proof of the second inequahty in (1671) let us recall that for any subsystem 
of d elements that we denote by Mi, ■ ■ ■ , the length perpendicular h from the origin to 
the affine hyperspace generated by the end points of vectors {u\, ■ ■ ■ , u^} is expressed, by 
dH, as 

h = ^^"^'---'^^^ ■ (68) 
\\R{ui, ■ ■ ■ ,Ud)\\ 

From the definition of D as a determinant it follows that 

d 

D{{zou,}) = l[ziD{{u,}). 

1=1 

Next, for the jth co-ordinate of the rotor R{{z o m}) we have 

1 

R'{{zou,}) = -l[ziR\{ui}), 



Zj 

^ 1=1 



so that 



\2 



n^^— < ^f[ziVd6{z)\\Ri{uM- 



mm-i Zj \z, 
1=1 ^ ^ II 1=1 



Hence 



\\R{{z o Ui})\\ 



implying the second inequality in (1^ . □ 
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6 Back to options; properties of solutions: non-expansion 
and homogeneity 

Let us now calculate the reduced Bellman operator of European colored options given by 
(fTOj) . Changing variables ^ = {C,^, . . . , ^"^) to r] = o z yields 

J 

{Bf){z\...,z'')= mm max [/(ry) - J] y (ry^ - pz^)], (69) 

1=1 



or, by shifting, 



(S/)(^\...,z^)=min ^ max Jf{v)-h,v)] (70) 

7 {?yG[2H<ii-p).-2'K-p)]} 



with /(?7) = f{rj + pz). Assuming / is convex (possibly not strictly), we find ourselves in 
the setting of Section [3] with 11 being the rectangular parallelepiped 

n.,p = y<i=i[z\di - p), z\ui - p)], 

with vertices 

Vi = ^1° z- pz, 

where 

= {di\iei, Ujlj^}, 
are the vertices of the normalized parallelepiped 

n= xtjd„M,], (71) 

parametrized by all subsets (including the empty one) / C {1, . . . , J}. 

Since the origin is an internal point of 11 (because di < p < Ui), condition (ii) of 
Theorem 14.11 is satisfied. Condition (i) is rough in the sense that it is fulfilled for an 
open dense subset of pairs {di,Ui). Applying Theorem 14.11 (and Remark [5]) to (170|) and 
returning back to ^ yields the following. 

Theorem 6.1. // the vertices C,i of the parallelepiped U are in general position in the 
sense that for any J subsets /i, ■ ■ ■ , Ij, the vectors {C,i^ — pl}'l^i are independent in , 
then 

iBf){z) = maxEnf{^oz), z = {z\- ■ ■ ,z'), (72) 

where {VL} is the collection of all subsets Q = C^j^, ■ ■ ■ of the set of vertices ofU, of 

size J + 1, such that their convex hull contains pi as an interior point (1 is the vector 
with all coordinates 1), and where denotes the expectation with respect to the unique 
probability law {pi}, G Vt, on the set of vertices o/ 11, which is supported on f2 and is 
risk neutral with respect to pi, that is 

/C{1,...,J} 

Moreover, if 

fii ° ^) - (7/„... a - Pl) o z) > fiC o z) - (77,,.. ,7,^,, (C - pi) o z) 
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for all vertices such that ^ G and ( ^ Q, where 'Jh,--- is the corresponding 
optimal value for the polyhedron n[^/j, ■ ■ ■ , ^/j+J, then 

{Bf){z\...,z'')=Enf{^oz). (74) 

Risk neutrality now corresponds to its usual meaning in finances, i.e. ( 1731) means that 
all discounted stock prices are martingales. 

Notice that the max in (172]) is over a finite number of explicit expressions, which is 
of course a great achievement as compared with initial minimax over an infinite set. In 
particular, it reduces the calculation of the iterations to the calculation for a control 
Markov chain. Let us also stress that the number of eligible fl in (I72|) is the number of 
different pyramids (convex polyhedrons with J + 1 vertices) with vertices taken from the 
vertices of 11 and containing pi as an interior point. Hence this number can be effectively 
calculated. 

Remark 9. Here we used the model of jumps, where each can jump independently in 
its interval. Thus we used the theory of Section\^ only for the case of a polyhedral H being 
a parallelepiped. The results of Section are given in a more general form to allow more 
general models of correlated jumps, see end of Section UU 

Let us point our some properties of the operator B given by ( 172]) that are obvious, but 
important for practical calculations: it is non-expansive: 

||S(/i)-i3(/2)||<||/l-/2||, 

and homogeneous (both with respect to addition and multiplication): 

B{X + f) = X + B{f), B{Xf) = XB{f) 

for a function / and A G R (resp. A > 0) for the first (resp second) equation. Finally, if 
fp is a power function, that is 

fp{z) = {zj^---{z'y^, 

then fp{$, oz) = fpiOfpiz) implying 

{B-fp)iz) = {{Bf,){l)rfp{z). (75) 

Therefore, power functions are invariant under B (up to a multiplication by a constant). 
Consequently, if for a payoff / one can find a reasonable approximation by a power 
function, that is there exists a power function fp such that ||/ — /p|| < e, then 

Wf-X^fp\\ = \\f-fp\\<e, X = {Bfp){l), (76) 

so that an approximate calculation of B"'f is reduced to the calculation of one number A. 
This implies the following scheme for an approximate evaluation of B: first find the best 
fit to / in terms of functions a + fp (where fp is a power function and a a constant) and 
then use (176]) . 
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7 Sub- modular payoffs 



One can get essential reduction in the combinatorics of Theorem 16.11 (i.e. in the number 
of ehgible Q) under additional assumptions on the payoff /. The most natural one in the 
context of options turns out to be the notion of sub-modularity. A function / : — > R+ 
is called sub-modular, if the inequality 

fixi, 1/2) + f{x2, yi) > fixi, xi) + f{yi, 7/2) 

holds whenever xi < yi and X2 < y2- Let us call a function / : H'^ — )■ R+ sub-modular ii 
it is sub-modular with respect to any two variables. 



Remark 10. // / is twice continuously differentiable, then it is sub-modular if and only 

dzidzi 



^/7^<0/or allt^j. 



As one easily sees, the payoffs of the first three examples of rainbow options, given in 
Section [2], that is those defined by ([1]), ([2]), (|3]), are sub-modular. Let us explain, on the 
examples of two and three colors J, how the assumptions of sub-modularity can simplify 
Theorem 16. 1[ 

Let first J = 2. The polyhedron flTTl) is then a rectangle. From sub- modularity of / it 
follows that if Q in Theorem 16.11 is either 

^12 = {{di, d2), (di, U2), (ui, U2)}, 

or 

^21 = {{di, d2), d2), U2)}, 

then (/, — (70,0 coincide for all vertices ^ of 11. Hence Q12 and Q21 can be discarded 
in Theorem 16.11 i.e the maximum is always achieved either on 

= {{di, d2), (rfi, U2), {ui, ^2)}, 

or on 

= {(c/i, M2), ^2), M2)}- 

But the interiors of the triangle formed by Qu and Qd do not intersect, so that each point 
of n (in general position) lies only in one of them (and this position does not depend any 
more on /). Hence, depending on the position of pi in H, the expression (!72|) reduces 
either to Ej^^ or to E^^. This yields the following result (obtained in pS]). 

Theorem 7.1. Let J = 2 and f be convex sub-modular. Denote 

^ _ {uiU2 - did2) - p{ui -di + U2- d2) p-di _ p - ^2 ^^^^ 

{ui - di){u2 - d2) ui-di U2-d2 

If K>0, then {Bf){zi,Z2) equals 

— — '^f{uizi, d2Z2) + - — ^-f{dizi, U2Z2) + Kf{diZi, d2Z2), (78) 
Ui — di U2 — "2 

and the corresponding optimal strategies are 

1 _ f{UiZi,d2Z2) - f{diZi,d2Z2) 2 _ f{dlZl,U2Z2) - f{diZi,d2Z2) 

^ zi{ui-di) ' ^ Z2{u2-d2) 
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If K <0, the {Bf){zi,Z2) equals 



^ f{diZi,U2Z2) + — ^f{uiZi,d2Z2) + \K\f{uiZi,U2Z2), (79) 



Ml - di U2- d2 

and 

1 _ f{UiZi,U2Z2) - f{diZi,U2Z2) 2 _ /(Wl^l,^2^2) - /(^l^l. ^^2^2) 

zi{ui - di) ' ^ Z2{u2-d2) 

Clearly the linear operator B preserves the set of convex sub-modular functions. Hence 
one can use this formula recursively to obtain all powers of in a closed form. For 
instance in case k = one obtains for the hedge price the following two-color extension 
of the classical Cox-Ross-Rubinstein formula: 

B-fiSlS^,) = p-J2c'n [{^J (^^) Mdr'Sld'^ur'S',). (80) 

Now let J = 3. Then polyhedron ( |7T1) is a parallelepiped in R^. From sub-modularity 
projected on the first two co-ordinates we conclude that whenever the vertices (c?i, d2, (is) 
and (mi, U2, d^) are in Q, then 

/(0-(7o,0 

should coincide for being (di, 6/2, rfs), {ui,U2,d3), (mi, 6/2, (is), (t/i, ^2, (^3). In other 
word, the pair ((ii, (^2, (is), (mi, ^2, (is) can be always substituted by the pair (-ui, c?2, (is), 
((ii, ^2, (is)- Consequently, those fl containing the pair ((ii, (i2, (is), (mi, M2, (is) are super- 
fluous, they can be discarded from the possible fl competing in formula (!72|) . Similarly, 
we can discard all those Q containing six pairs, three of which containing {di, c?2, (is) and 
one among ((ii, (i2, Ms); (c^i, "^2, (is), (mi, (i2, c^s), and other three containing (mi,M2,M3) and 
one among ((ii,U2,Ms), (ui,U2,c?3), (ui, (i2, ^^s)- 

These considerations reduce dramatically the number of eligible Q. In particular, if 
pi lies in the tetrahedron Qd formed by the vertices ((ii, (^2, (is), ((ii, (i2, its), ((ii, ^i2, c?s), 
(mi, (i2, (is), then the only eligible Q is fid- If pl lies in the tetrahedron flu formed by 
the vertices (mi, U2, U3), {di, U2, us), (mi, U2, (^3), {ui, d2, Ms), then the only eligible fl is flu- 
Formally these cases are easily seen to be distinguished by the inequalities ci;i23 > and 
tti2s < — 1 respectively, where 

/ Ml - p U2- p u^- p 

"123 = 1 - 



Ml — di M2 — C?2 Ms — C?3 

This yields the following result (by continuity we are able to write ^123 > instead of 
a strict inequality), where we use the following notation: for a set / C {1, 2, J}, fi{z) 
is fi^^Zi, ■ ■ ■ , ^"^Zj) with = di for i E I and = m, for i ^ I. 

Theorem 7.2. Let J = 3 anc? / be continuous convex and sub-modular, 
(i) If ai23 > 0, then 



m{z) 



"123/0 (^) + 



Ml - p 

Ml — di 



/{1}(^) + 



"^"'/{2}(^) + :^/{3>(.) . (82) 



M2 - d2 



Ms - C?3 



(a) If ai23 < —1, then 
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(83) 



Hence in these cases, our B again reduces to a linear form, allowing for a straightfor- 
ward calculation of its iterations, as in case J = 2 above. 

Suppose now that pi lies neither in the tetrahedron VL^., nor in VL^ (i.e. neither of 
the conditions of Theorem 17.21 are satisfied). From the above reductions of possible f2, it 
follows that in that case one can discard all Vl containing either (rfi, ^2, d^) or (-ui, M2, ms). 
Hence only six vertices are left for eligible VL. From the consideration of general position 
we further deduce that altogether only six Vt are possible, namely the three tetrahedrons 
containing the vertices (rfi, (i2, M3), (lii, M2, c^s), (mi, (^2, dz) and one vertex from (c^i, M2, M3), 
(mi, -^2, (is), (Mi,(i2,M3), and symmetrically the three tetrahedrons containing the vertices 
(c?i, M2, Ms), (mi, M2, c^s), («!, c?2, ^s) and one vertex from (rfi, ^2, Ms), {di-, M2, c^s), (^i, c?2, d^). 
However, any particular point in general position belongs to only three out of these six 
leaving in formula 0721) the max over three possibilities only. The particular choice of 
these three tetrahedrons depends on the coefficients 




u\—r U2—r 



Ul- 


-di 


U2~ 


~d2 


"1 




"3 


—r 


Ml- 


-di 


M3- 


-ds 


U2- 


-r 


"3- 


-r 



U2—d2 Uz—d3 



and leads to the following result obtained in Hucki and Kolokoltsov |20] (though with 
much more elaborate proof than here). 

Theorem 7.3. Let again f be convex and sub-modular, but now > 0123 > — 1- 

(i) If a 12 > 0, > and 0:23 > 0, then 

( (-«123) /{1,2}(Z) + ai3/{2}(z) + a23/{l}(z) + ^^/{3}(z) 

(i3/)(z) = i max <^ (-^123) /{i,3}(z) + ai2/{3}(z) + a23/{i}(z) + ^J{2}{z) 

[ (-"123) /{2,3}(z) + Q;i2/{3}(z) + Q;13/{2}(z) + ^^/{l}(z) 

(ii) If aij < 0, ttjk > and > 0, 

where {i,j,k} is an arbitrary permutation of the set {1,2,3}, then 

{(-tty-fc) /{i,i}(z) + «ifc/{i}(z) + ajkf{i}{z) + ^J^/{fc}(z) 
",fe/{i}(z) + (-a.,)/{.,,}(z) + ^J{^Mi^) - 
a^fe/b■}(z) + (-ai,-)/{i,i}(z) + ^/{i,fc}(z) - ^/w(z) 
(Hi) If aij > 0,ajk < and a^k < 0, 

where {i,j, k} is an arbitrary permutation of the set {1, 2, 3}, then 

( «u/{fc}(z) + (-a,fc)/{,,fe}(z) + ^J{i,k}i^) - S/b}(z) 

(i3/)(z) = i max <^ aijf{k}{z) + (-a,fe)/{i,fc}(z) + ^f{j,k}i^) - ^/«(z) 
[ ("123 + l)/{fc}(z) - ajkf{j,k}iz) - aikf{i,k}iz) - ^Jfi^/{i,j}(z) 

One has to stress here that the application of Theorem 17.31 is rather limited: as B is 
not reduced to a linear form, it is not clear how to use it for the iterations of B, because 
the sub-modularity does not seem to be preserved under such B. 
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8 Transaction costs 

Let us now extend the model of Section [2] to include possible transaction costs. They can 
depend on transactions in various way. The simplest for the analysis are the so called 
fixed transaction costs that equal to a fixed fraction (1 — /3) (with /3 a small constant) of 
the entire portfolio. Hence for fixed costs, equation ([6]) changes to 

J J 

As one easily sees, including fixed costs can be dealt with by re-scaling p, thus bringing 
nothing new to the analysis. 

In more advanced models, transaction costs depend on the amount of transactions 
(bought and sold stocks) in each moment of time, i.e. are given by some function 

and are payed at time, when the investor changes 7m-i to 7^- In particular, the basic 
example present the so called proportional transaction costs, where 

J 
i=i 

(again with a fixed /3 > 0). We shall assume only that g has the following Lipshitz 
property: 

\9hi,z) - 5'(72,2;)| < I3\z\\-fi - 72I (86) 

with a fixed (3 > 0. 

To deal with transaction costs, it is convenient to extend the state space of our game, 
considering the states that are characterized, at time m — 1, by 2 J + 1 numbers 

When, at time m — 1, the investor chooses his new control parameters 7^, the new state 
at time m becomes 

where the value of the portfolio is 

J J 

X„, = ^ lln^lnSln-l + P(^m-1 " ^ llnSL-l) " dilm " ^m-l, ^m-l)- (87) 

The corresponding reduced Bellman operator from Section [2] takes the form 

{Bf){z,v) = minmax[/(^o2;,7) - (^^^oz-pz) +g{-f-v,z)], (88) 

7 ? 

where z, f G R'^, or, changing variables ^ = {^^, . . . , ^'^) to rj = ^ o z and shifting, 

(^/)(^,^) = min max [f{r] + pz,-^)-{-^,ri) + g{-f-v,z)]. (89) 

7 WelzJ{d,-p),zi{u,-p)]} 
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On the last step, the function / does not depend on 7, so that Theorem 15.11 can be used 
for the calculation. But for the next steps Theorem 15.21 is required. 
For its recursive use, let us assume that 

\fiz,vi) - fiz,V2)\ < a\z\\vi -V2\, 

and a is small enough so that the requirements of Theorem 15.21 are satisfied for the r.h.s. 
of OHS]). By Theorem E21 

iBf){z, v) = maxEf,[/(e o z, 7^) +g{ln- v, z)]. (90) 

Notice that since the term with v enters additively, they cancel from the equations for 
7n, so that the values of 7q do not depend on v. Consequently, 

\iBf)iz,v,) - iBf){z,V2)\ < m^xEn[gi7n-Vi,z)-g{^n-V2),z)] < (3\z\\v, - V2\. (91) 

Hence, if at all steps the application of Theorem 15.21 is allowed, then {B^ f){z,v) 
remains Lipshitz in v with the Lipshitz constant I3\z\ (the last step function does not 
depend on v and hence trivially satisfies this condition). 

Let Xi, X2 be the characteristics, defined before Theorem 15. 2[ of the set of vertices 
^/ — pi of the parallelepiped ^'j=i[dj,Uj\ — pi. By Proposition 15. H the corresponding 
characteristics xi(z), ^2(2;) of the set of vertices of the scaled parallelepiped 

yij^i[z^dj, z^Uj] — pz 

have the lower bounds ^ 

Xi{z) > \z\Xi——, 2 = 1,2. 

dd{z) 

As in each step of our process the coordinates of z are multiplied by dj or Uj, the corre- 
sponding maximum 6n{z) of the 6 of all z that can occur in the n-step process equals 

5„(z) = 5(.)f^^V. (92) 



Thus we arrive at the following result. 
Theorem 8.1. Suppose (3 from (186|1 satisfies the estimate 

1 



(3 < min(xi, X2) 



d5Jz) 



where 5n{,z) is given by (p2|) . Then the hedge price of a derivative security specified by 
a final payoff f and with transaction costs specified above is given by ( !T2|) . where B is 
given by ( 188|) . Moreover, at each step, B can be evaluated by Theorem \5.^ i.e. by ( 190|) . 
reducing the calculations to finding a maximum over a finite set. 

Of course, for larger /3, further adjustments of Theorem 15.21 are required. 
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9 Rainbow American options and real options 



In the world of American options, when an option can be exercised at any time, the 
operator BG'(X, S^, S"^) from ([9]) changes to 

BG{X,S\--- ,3-') 



max mm 

7 



J J 

G{X, S\--- ,S'),-min G{pX + ^ 7^^' - P l'S\ ■ ■ ■ , ^S') 



i=l i=l 



(93) 



so that the corresponding reduced operator takes the form 
{Bf){z^, z"^) = minmax 



7 



i=l 



(94) 



or equivalently 

{Bf){z^, z'^) = max 



p/(pz),minmax[/(e^;.\eV, ■ ■ ■ ,eV) - J^f ^^(e - p)] 
7 5 



i=l 



Consequently, in this case the main formula (I72|) of Theorem 16.11 becomes 



z'^) = max 



pf{pz),maxEQf{^o z) 



(95) 



(96) 



which is of course not an essential increase in complexity. The hedge price for the n-step 
model is again given by (fT2|) . 

Similar problems arise in the study of real options. We refer to Dixit and Pindyck 
[16] for a general background and to Bensoussan et al [7] for more recent mathematical 
results. A typical real option problem can be formulated as follows. Given J instruments 
(commodities, assets, etc), the value of the investment in some project at time m is 
supposed to be given by certain functions fm{S^,--- , S^) depending on the prices of 
these instruments at time m. The problem is to evaluate the price (at the initial time 0) 
of the option to invest in this project that can be exercised at any time during a given 
time-interval [0,T]. Such a price is important, since to keep the option open a firm needs 
to pay ceratin costs (say, keep ready required facilities or invest in research). We have 
formulated the problem in a way that makes it an example of the general evaluation of 
an American rainbow option, discussed above, at least when underlying instruments are 
tradable on a market. For practical implementation, one only has to keep in mind that the 
risk free rates appropriate for the evaluation of real options are usually not the available 
bank accounts used in the analysis of financial options, but rather the growth rates of 
the corresponding branch of industry. These rates are usually estimated via the CAPM 
(capital asset pricing model), see again [16]. 



10 Path dependence and other modifications 

The Theory of Section E] is rough, in the sense that it can be easily modified to accom- 
modate various additional mechanisms of price generations. We have already considered 
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transaction costs and American options. Here we shall discuss other three modifications: 
path dependent payoffs, time depending jumps (including variable volatility) and non- 
linear jump formations. For simplicity, we shall discuss these extensions separately, but 
any their combinations (including transaction costs and American versions) can be easily 
dealt with. 

Let us start with path dependent payoffs. That is, we generalize the setting of Section 
[2] by making the payoff / at time m to depend on the whole history of the price evolutions, 
i.e. being defined by a function f{So, ■ ■ ■ , Sm), Si = {S^ ■ ■ ■ , 5/), on R-^('"+i). The 
state of the game at time m must be now specified by (m + 1) J + 1 numbers 

X„, Si = {Si, ■ ■ ■ , 5/), i = 0, ■ ■ ■ , m. 

The final payoff in the n-step game is now G = X — f{So, ■ ■ ■ , Sn) and at the pre ultimate 
period n — 1 (when 5*0, ■ ■ ■ , 5^-1 are known) payoff equals 

BG(X, 5*0, ■ ■ ■ ,5'„_i) = X - -minmax[/(5'o, ■ ■ ■ ,Sn-i,C° S^-i) - (7, '5'„_i o (^ - pi))] 

pi £. 

= X (S„_i/)(S'o, ■ ■ ■ , 5'„„i), 

P 

where the modified reduced Bellman operators are now defined as 
(i3^_i/)(zo, ■■■ ,^m-i) = min max [f [zq,- ■■ ,2^-1,^0 Zm-i) - [l.^o z - pz)]. (97) 

Consequently, by dynamic programming, the guaranteed payoff at the initial moment 
of time equals 

X - 1^0(^1 ■■■(^„-i/)---), 

and hence the hedging price becomes 

i7" = 4^i3o(i5i---(i3„_i/)---). (98) 

No essential changes are required if possible sizes of jumps are time dependent. Only 
the operators Bm-i from fl971) have to be generalized to 

(i3„^_i/)(2:o, ■ ■ ■ , Zm_i) = min max [f{zQ,---,Zm-i,ioZm-i)-{l,ioz-pz)], (99) 

where the pairs (rf™', m™), j = 1, ■ ■ ■ , J, m = 1, ■ ■ ■ , n specify the model. 

Let us turn to nonlinear jump patterns. Generalizing the setting of Section [2] let us 
assume, instead of the stock price changing model S^+i = ° Sm, that we are given k 
transformations Qi : R"' — t- R"', z = 1, ■ ■ ■ , fc, which give rise naturally to two models of 
price dynamics: either 

(i) at time m + 1 the price S^+i belongs to the closure of the convex hull of the set 
{5'i(5'm)}, i = 1, ■ ■ ■ ,k (interval model), or 

(ii) Sm+i is one of the points {gi{Sm)}, 2 = 1, ■ ■ ■ , fc. 

Since the first model can be approximated by the second one (by possibly increasing 
the number of transformations Qi), we shall work with the second model. 
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Remark 11. Notice that maximizing a function over a convex polyhedron is equivalent 
to its maximization over the edges of this polyhedron. Hence, for convex payoffs the two 
models above are fully equivalent. However, on the one hand, not all reasonable payoffs 
are convex, and on the other hand, when it comes to minimization (which one needs, say, 
for lower prices, see Section [77]) . the situation becomes rather different. 

Assuming for simplicity that possible jump sizes are time independent and the payoffs 
depend only on the end-value of a path, the reduced Bellman operator (fTU|) becomes 

{Bf){z)=mm max [f{g,{z))-{^,g,{z)~pz)], z = {z\ z'), (100) 

7 iG{l, ■■■,*:} 

or equivalently 

{Bf){z) = min max [f{r]i + pz) - (7, 77^)]. (101) 

The hedge price is still given by f lT2|) and operator f llOOp is calculated by Theorem 13.21 
It is worth noting that iik = d+1 and {gi{z)} form a collection of vectors in a general 
position, the corresponding risk-neutral probability is unique. Consequently our hedge 
price becomes fair in the sense of 'no arbitrage' in the strongest sense: no positive surplus 
is possible for all paths of the stock price evolutions (if the hedge strategy is followed). 
In particular, the evaluation of hedge strategies can be carried out in the framework of 
the standard approach to option pricing. Namely, choosing as an initial (real world) 
probability on jumps an arbitrary measure with a full support, one concludes that there 
exists a unique risk neutral equivalent martingale measure, explicitly defined via formula 
( l33l) . and the hedge price calculated by the iterations of operator f llOip coincides with the 
standard risk-neutral evaluation of derivative prices in complete markets. 

11 Upper and lower values; intrinsic risk 

The celebrated non-arbitrage property of the hedge price of an option in CRR or Black- 
Scholes models means that a.s., with respect to the initial probability distribution on 
paths, the investor cannot get an additional surplus when adhering to the hedge strategy 
that guarantees that there could be no loss. In our setting, even though our formula in 
case J = 1 coincides with the CRR formula, we do not assume any initial law on paths, 
so that the notion of 'no arbitrage' is not specified either. 

It is a characteristic feature of our models that picking up an a priori probability 
law with support on all paths leads to an incomplete market, that is to the existence 
of infinitely many equivalent martingale measures, which fail to identify a fair price in 
a unique consistent way. Notice that our extreme points are absolutely continuous, but 
usually not equivalent to a measure with a full support. 

Remark 12. The only cases of a complete market among the models discussed above are 
those mentioned at the end of Section [721 that is the models with precisely J + 1 eligible 
jumps of a stock price vector S in each period. 

For the analysis of incomplete markets is of course natural to look for some subjective 
criteria to specify a price. Lots of work of different authors were devoted to this endeavor. 
Our approach is to search for objective bounds (price intervals), which are given by our 
hedging strategies. 
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Remark 13. Apart from supplying the lower price (as below), one can also argue about 
the reasonability of our main hedging price noting that a chance for a possible surplus can 
be (and actually is indeed) compensated by inevitable inaccuracy of a model, as well as 
by transaction costs (if they are not taken into account properly). Moreover, this price 
satisfies the so called 'no strictly acceptable opportunities' (NSAO) condition suggested in 
Carr, Geman and Madan J7^ . 

For completeness, let us recall the general definitions of lower and upper prices, in 
the game theoretic approach to probability, given in Shafer and Vovk [SU]. Assume a 
process (a sequence of real numbers of a fixed length, say n, specifying the evolution of 
the capital of an investor) is specified by alternating moves of two players, an investor 
and the Nature, with complete information (all eligible moves of each player and their 
results are known to each player at any time, and the moves become publicly known at the 
moment when decision is made). Let us denote by X°(^) the last number of the resulting 
sequence, starting with an initial value a and obtained by applying the strategy 7 (of the 
investor) and ^ (of the Nature). By a random variable we mean just a function on the set 
of all possible paths. The upper value (or the upper expectation) E/ of a random variable 
/ is defined as the minimal capital of the investor such that he/she has a strategy that 
guarantees that at the final moment of time, his capital is enough to buy /, i.e. 

E/ = inf{« : 37 : V^, X-(0 - /(O > 0}. 

Dually, the lower value (or the lower expectation) E/ of a random variable / is defined 
as the maximum capital of the investor such that he/she has a strategy that guarantees 
that at the final moment of time, his capital is enough to sell /, i.e. 

E/ = sup{« : 37 : V^, X^"(0 + /(O > 0}. 

One says that the prices are consistent if E/ > E/. If these prices coincide, we are in a 
kind of abstract analog of a complete market. In the general case, upper and lower prices 
are also referred to as a seller and buyer prices respectively. 

It is seen now that in this terminology our hedging price for a derivative security is the 
upper (or seller) price. The lower price can be defined similarly. Namely, in the setting 
of Section |51 lower price is given by 

^(^Io«;/)('S'd, ■ ■ ■ , 5*0 ), 

where 

{Biou,f){z) =m&i^ niin [/(^ o 2;) - (7, ^ o 2; - p^;)]. (102) 

In this simple interval model and for convex / this expression is trivial, it equals f{pz). 
On the other hand, if our / is concave, or, more generally, if we allow only finitely many 
jumps, which leads, instead of fll02p . to the operator 

(i3«o«;/)(^)=max min [/(^ o 2;) - (7, ^ o 2; - p^;)], (103) 

then Theorem 14.31 applies giving for the lower price the dual expression to the upper price 
(1721) . where maximum is turned to minimum (over the same set of extreme risk- neutral 
measures) : 

{Bio^f){z)=mmY.nf{ioz), z = {z\- ■ ■ ,z'). (104) 



31 



The difference between lower and upper prices can be considered as a measure of 
intrinsic risk of an incomplete market. 



12 Continuous time limit 

Our models and results are most naturally adapted to discrete time setting, which is not 
a disadvantage from the practical point of view, as all concrete calculations are anyway 
carried out on discrete data. However, for qualitative analysis, it is desirable to be able to 
see what is going on in continuous time limit. This limit can also be simpler sometimes 
and hence be used as an approximation to a less tractable discrete model. Having this 
in mind, let us analyze possible limits as the time between jumps and their sizes tend to 
zero. 

Let us work with the general model of nonlinear jumps from Section [TOl with the re- 
duced Bellman operator of form f llOOp . Suppose the maturity time is T. Let us decompose 
the planning time [0,T] into n small intervals of length t = T/n, and assume 

= ^ + r-0,(z), ^ = l,---,fc, (105) 

with some functions 0j and a constant a G [1/2, 1]. Thus the jumps during time r are of 
the order of magnitude r". As usual, we assume that the risk free interest rate per time 
r equals 

p = 1 + rr, 

with r > 0. 

From ( llOOp we deduce for the one-period Bellman operator the expression 

Brf{z) = ^^max Vp[(z,r)/(z + r>,(^)), (106) 

where / are subsets of {1, ■ ■ ■ ,n} of size |/| = J + 1 such that the family of vectors 
z + T^cpilz), i E I, are in general position and {pI{z,t)} is the risk neutral probability 
law on such family, with respect to pz, i.e. 

J2pi {z, t){z + r>,(^)) = (1 + rT)z. (107) 

Let us deduce the HJB equation for the limit, as r — )■ 0, of the approximate cost- 
function B^-\ t e [0, T], with a given final cost /t, using the standard (heuristic) dynamic 
programming approach. Namely, from (11061) and assuming appropriate smoothness of / 
we obtain the approximate equation 



i&I 



Mz) + r-^Uz) + Ir^- i^Hz), ] + 0{t^-' 



dz 2 \ dz"^ 

Since {pi } are probabilities and using (I107P , this rewrites as 
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where 

r— )-0 

(clearly well defined non-negative numbers). This leads to the following equations: 



in case a = 1/2, and to the trivial first order equation 

df ^ . df 



rf = % + r{zM) (109) 



with the obvious solution 

/(t,z) = e-^(^-*)/r(e-^(^-*)^) (110) 

in case a > 1/2. 

Equation (11081) is a nonlinear extension of the classical Black-Scholes equation. Well 
posedness of the Cauchy problem for such a nonlinear parabolic equation in the class of 
viscosity solutions is well known in the theory of controlled diffusions, as well as the fact 
that the solutions solve the corresponding optimal control problem, see e.g. Fleming and 
Soner [17] . 

Remark 14. Having this well posedness, it should not be difficult to prove the conver- 
gence of the above approximations rigorously, but I did not find the precise reference. 
Moreover, one can be also interested in path-wise approximations. For this purpose a 
multidimensional extension of the approach from Bick and Willinger fW^ (establishing 
path-wise convergence of Cox-Ross-Rubinstein binomial approximations to the trajecto- 
ries underlying the standard Black-Scholes equation in a non-probabilistic way) would be 
quite relevant. 

In case J = 1 and the classical CCR (binomial) setting with 

a/t^i = {u — 1)z = o-Jtz, \fr<^2 = (d — l)z = —o\frz, 

equation (11081) turns to the usual Black-Scholes. 

More generally, if /c = J + 1, the corresponding market in discrete time becomes 
complete (as noted at the end of Section [10]) . In this case equation (I108P reduces to 

which is a generalized Black-Scholes equation describing a complete market (with ran- 
domness coming from J correlated Brownian motions), whenever the diffusion matrix 

J+i 



{P'),^ = Y.V^{ZW^^Z),<^\{Z). 
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As a more nontrivial example, let us consider the case of J = 2 and a sub-modular final 
payoff /t, so that Theorem 17.11 applies to the approximations Br- Assume the simplest 
(and usual) symmetric form for upper and lower jumps (further terms in Taylor expansion 
are irrelevant for the limiting equation): 

Ui = l + ai^fr, di = \- Oiy/r, i = 1,2. (112) 

Hence 

Ui-p 1 r . 

/r, 1 = 1,2, 



Ui — di 2 2ai 
and 

1^/1 In 
= -7;rVT{— + —). 

2 (Jl (72 

As K < 0, we find ourselves in the second case of Theorem 17.11 Hence the only eligible 
collection of three vectors is (di, ti2), (mi, ^2), (wi, ^2), and the probability law p[ is 
(1/2, 1/2,0). Therefore, equation f llOSp takes the form 



r 9f , df, 1 



2 2d'f . d'f , 2 2d'f 



(113) 



The limiting Black-Scholes type equation is again linear in this example, but with 
degenerate second order part. In the analogous stochastic setting, this degeneracy would 
mean that only one Brownian motion is governing the behavior of both underlying stocks. 
This is not surprising in our approach, where Nature was assumed to be a single player. 
One could expect uncoupled second derivatives (non-degenerate diffusion) in the hmit, if 
one would choose two independent players for the Nature, each playing for each stock. 

Thus we are still in the setting of an incomplete market. The hedge price calculated 
from equation flll3p is actually the upper price, in the terminology of Section [TTJ To get 
a lower price, we shall use approximations of type ( I104p . leading, instead of (llOSp . to the 
equation 

If J = 2 the payoff is submodule, the maximum can be taken over the triples {di, ^2), 
{di,U2), {ui,d2) or (lii, M2), (mi, ^2), (mi, "^2) (under flll2p only the second triple works). 
Similarly the minimum can be taken only over the triples (cii, (^2), ((ii, M2), (^^i, ^2) or 
((^1,^2), {ui,d2), {ui,U2). Under flll2p both these cases give the same limit as r — )■ 0, 
yielding for the lower price the equation 



X df , df, 1 
rf = — — h r{z, ——) H — 
^ dt ^ ' dz' 2 



2 25V , 2 25V 

+ 2^,^2.1.2^^ + ^2^2^ 



2 J 



(115) 



that differs only by sign at the mixed derivative from the equation for the upper price. 

As / was assumed sub-modular, so that its mixed second derivative is negative, we 
have 

(Tl(T2^ 7^ < < -(7lO-2T^ 

OZ1OZ2 OZ1OZ2 

Hence, for the solution fu of the upper value equation f lll3p . the solution fi of the lower 
value equation f lllSp . and the solution /c of the classical Black-Scholes equation of a 
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complete market based on two independent Brownian motions, i.e. equation flllSp or 
(11151) without the term with the mixed derivative (with the same sub-modular initial 
condition fx) we have the inequality 

fl<fc< fu, 

as expected. 

Equations (11131) and (I115p can be solved explicitly via Fourier transform, just as the 
standard Black-Scholes equation. Namely, changing the unknown function f to g hj 



so that 



fizi,Z2 

df 



e '■^^ ^^g{—\ogzi, — \ogZ2) 
Ox 02 



e ^(^-) — l^(-log.i,Alog.,), 



dzi CTiZidyi ai 

transforms these equations to the equations 



0-2 



d'^q d'^q d'^q 

± ip 2 - I — 

dyl dyidy2 dyl 







(116) 



(with =F respectively). Equation (I116p has constant coefficients and the equation for the 
Fourier transform g{p) of g is obviously 



= -j[{PiTP2f - iC^r - ai)pi - i{2r - ^2)^2]^- 



:il7) 



Hence the inverse Cauchy problem for equation (11160 with a given final function gT equals 
the convolution of gr with the inverse Fourier transform of the functions 



exp{--(T - t)[{pi TP2Y - ^(2r - ai)pi - i{2r - 0-2)^2]}, 
which equal (after changing the integration variables pi and p2 to qi = pi—p2, q2 = P1+P2) 

/ rfgirfg2exp{-^(T-t)g2 



2(27 

(ai -a2)(T-t) 



+ '-f{yi + y2 + {2r 



ai + ^2 



){T-t)]} 



(with gi_2 corresponding to =f), or exphcitly 



1 



2 v/27r(r - t) 
and 



yi + y2 , , 01 + 02. , r 



KT - 1) 



yi - y2 



(01 -02)(T-t) 



1 



1 



2 ^2n{T - t) 



6 



yi - y2 _ (oi - 02)(r - 1) 
2 4 



exp{- 



KT-t) 



l/i + l/2 + (2r-^i±^)(r-t) 
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respectively, where 6 denotes the Dirac 5-function. Returning to equations f lllSp and 
(11151) we conclude that the the solutions /„ and // respectively of the inverse time Cauchy 
problem for these equations are given by the formula 

/»00 /"OO 

fu,lit,Zi,Z2) = / G^_^{Zi,Z2;Wi,W2)fT{Wi,W2)dWidW2, (118) 

Jo Jo 

with the Green functions or transition probabilities being 

G:j,_^{zi,Z2;wi,W2) 

-5 - — log h - — log h (r )(T - t) 



and 



2n{T -t)(7ia2WiW2 K'^'^i W2 4 

— 7\ — log log o ^ 

8(T - 1) V(Ti wi a2 W2 2 J 

G^^^{zi,Z2;wi,W2) 

e-^(^-*) 1 ^1 1 , ^2 (ai-a2)(T-t) 
-d - — log — log 



2A/27r(T - t)(Tia2WiW2 K'^'^i 2a2 W2 4 

exp{--^ (- log ^ + 1 log ^ + (2r - ^)(T - t)) '} (120) 
8{T-t) \ai wi a2 W2 2 J 

respectively. Of course, formulas (IllSp can be further simplified by integrating over the 6- 
function. Singularity, presented by this 5-function, is due to the degeneracy of the second 
order part of the corresponding equations. 



13 Transaction costs in continuous time 

The difficulties with transaction costs are well known in the usual stochastic analysis 
approach, see e.g. Soner et al [52] and Bernard et al [llj. 

In our approach. Theorem 18.11 poses strong restrictions for incorporating transaction 
costs in a continuous limit. In particular, assuming jumps of size r° in a period of length 
r, i.e. assuming (11051) . only a = 1 can be used for the limit r — 0, because 5n{z) is of 
order (1 + r")", which tends tooo, asr = T/n— t-O, whenever a < 1. We know that for 
vanishing costs, assuming a = 1 leads to the trivial limiting equation (11091) . which was 
observed by many authors, see e.g. Bernard [9], McEneaney [H], Olsder [16]. However, 
with transaction costs included, the model with jumps of order r becomes not so obvious, 
but leads to a meaningful and manageable continuous time limit. To see this, assume 
that we are in the setting of Section [12] and transaction cost are specified, as in Section [8], 
by a function g satisfying (1861) . To write a manageable approximation, we shall apply the 
following trick: we shall count at time rm the transaction costs incurred at time r(m + 1) 
(the latter shift in transaction costs collection does not change, of course, the limiting 
process). Instead of (11061) we then get 
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[f{z + r>,(2;)) + g{j{z + r^z), r) - j{z, t),z + r0,(z))] , (121) 

where 7(z,r) is the optimal 7 chosen in the position z. Assuming g is different iable, 
expanding and keeping the main terms, yields the following extension of equation (11091) : 

r/ = ^ + r(z,^)+^(z), (122) 



where 



iG/ m,j=l 



with 7(2;) = \imr-^o-f{z,T). 

This is a non-homogeneous equation with the corresponding homogeneous equation 
being (11091) . Since the (time inverse) Cauchy problem for this homogeneous equation has 
the explicit solution (IllOp . we can write the explicit solution for the Cauchy problem of 
equation (11221) using the standard Duhamel principle (see e.g. [31]) yielding 

f{t,z) = e-''^^-'^fT{e-''^^-'^z) + e-'^'-'^tlj{e-'-^'-'h)ds. (123) 

The convergence of the approximations Sr^^'/r to this solution of equation (11221) follows 
from the known general properties of the solutions to the HJB equations, see e.g. [34] . 

Of course, one can also write down the modified equation (11081) obtained by introducing 
the transaction costs in the same way as above. It is the equation 

'^Z = ^ + ^) + 2 max (z) 



— ^,y.j,^,y.j J ^ ^ —^{-f{z))-^(j)i{z) 
^ m,j=l ' 



(124) 



However, as already mentioned, due to the restrictions of Theorem 18. H only the so- 
lutions to a finite difference approximations of equation (11241) (with bounded below time 
steps r) represent justified hedging prices. Therefore our model suggests natural bounds 
for time-periods between re-locations of capital, when transaction costs remain amenable 
and do not override, so-to-say, hedging strategies. Passing to the limit r — > in this model 
(i.e. considering continuous trading), does not lead to equation (I124p . but to the trivial 
strategy of keeping all the capital on the risk free bonds. This compelled triviality is of 
course well known in the usual stochastic setting, see e. g. Soner, Shreve and Cvitanic 



14 Fractional dynamics 

Till now we have analyzed the models where the jumps (from a given set) occur with 
regular frequency. However, it is natural to allow the periods between jumps to be more 
flexible. One can also have in mind an alternative picture of the model: instead of 
instantaneous jumps at fixed periods one can think about waiting times for the distance 
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from a previous price to reach certain levels. It is clear then that these periods do not 
have to be constant. In the absence of a detailed model, it is natural to take these waiting 
times as random variables. In the simplest model, they can be i.i.d. Their intensity 
represents a kind of stochastic volatility. Slowing down the waiting periods is, in some 
sense, equivalent to decreasing the average jump size per period. 

Assume now for simplicity that we are dealing with 2 colored options and sub-modular 
payoffs, so that Theorem 17. ll applies yielding a unique eligible risk- neutral measure. Hence 
the changes in prices (for investor choosing the optimal 7) follow the Markov chain X^{z) 
described by the recursive equation 

where (j){z) is one of three points {z^di, z'^U2), {z^ui, z'^d2), {z^ui, z'^U2) that are chosen 
with the corresponding risk neutral probabilities. As was shown above, this Markov chain 
converges, as r — > and n = [t/r] (where [s] denotes the integer part of a real number s), 
to the diffusion process Xf solving the Black-Scholes type (degenerate) equation (11131) . 
i.e. a sub-Markov process with the generator 



df 1 

Lfix) = -rf + r{z, -g^) + ^ 



(125) 



Assume now that the times between jumps Ti, T2, ■ ■ ■ are i.i.d. random variables with 
a power law decay, that is 

with P G (0, 1) (where P denotes probability and ~ means, as usual, that the ratio 
between the l.h.s. and the r.h.s. tends to one as t — ?■ 00). It is well known that such Tj 
belong to the domain of attraction of the /3-stable law (see e.g. Uchaikin and Zolotarev 
|55j ) meaning that the normalized sums 

e[ = r''(Ti + ■■■ + %.]) 

(where [s] denotes the integer part of a real number s) converge, as r — > 0, to a /3-stable 
Levy motion Qt, which is a Levy process on R+ with the fractional derivative of order /3 
as the generator: 

■^m - -^m - -fp^ f (/(* + r) - 

We are now interested in the process Y'^{z) obtained from X^{z) by changing the 
constant times between jumps by scaled random times Tj, so that 

Y;iz)=XJ,.iz), 

where 

iV; = max{u : < t}. 

The limiting process 

Nt = ma.x{u : 0^ < t} 

is therefore the inverse (or hitting time) process of the /3-stable Levy motion 0^. 

By Theorem 4.2 and 5.1 of Kolokoltsov [3^, (see also Chapter 8 in [31]), we obtain 
the following result. 
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Theorem 14.1. The process converges (in the sense of distribution on paths) to the 
process Yt = X^^, whose averages f(T — t,x) = Eif{YT-t{x)) , for continuous bounded f, 
have the explicit integral representation 

/»00 /»00 POO 

f{T-t,x)= / / / G~{zi, Z2]Wi,W2)Q(T - t,u) dudwidw2, 
Jo Jo Jo 

where , the transition probabilities of Xt, are defined by f lllOp . and where Q{t,u) 
denotes the probability density of the process Nt. 

Moreover, for f G C^(R'^), f{t,x) satisfy the (generalized) fractional evolution equa- 
tion (of Black-Scholes type) 

d^ t-^ 
-^/(«..) = L/(,x) + ^^M.). 

Remark 15. Similar result to Theorem 4-2 of Kolokoltsov [30] used above, but with 
position independent random walks, i.e. when L is the generator of a Levy process, were 
obtained in Meerschaert and Scheffler JJ^ , see also related results in Kolokoltsov, Korolev 
and Uchaikin JSMj, Henry, Langlands and Straka JJEl and in references therein. Rather 
general fractional Cauchy problems are discussed in Kochubei [27] . 

Similar procedure with a general nonlinear Black-Scholes type equation fllOSp will 
lead of course to its similar fractional extension. However, a rigorous analysis of the 
corresponding limiting procedure is beyond the scope of the present paper. 
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